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Abstract 

There are no upper limits on the possible number of massive, singlet (right-handed) 
neutrinos that may participate in the seesaw mechanism, and some string construc- 
tions motivate seesaw models with up to 0(100) right-handed neutrinos. In this case, 
the seesaw mass scale can be significantly higher than that in the traditional scheme 
with just 3 right-handed neutrinos. We consider the possible phenomenological impli- 
cations of such models, in particular, for lepton-flavour violation and electric dipole 
moments. Since the neutrino masses depend on the Majorana mass scale linearly, 
while supersymmetric loop corrections depend on it logarithmically, the magnitude 
of lepton-flavour- and CP-violating transitions may increase with the multiplicity of 
the right-handed neutrinos and may be enhanced by orders of magnitude. We also 
point out that, in the context of leptogensis, the bounds on the reheating tempera- 
ture and the lightest neutrino mass get relaxed compared to those in the case of 3 
right-handed neutrinos. 



1 Introduction 



The seesaw mechanism is arguably the most attractive way to explain the smallness 
of neutrino masses [1]- [4]. In its conventional form, the seesaw invokes 3 heavy sin- 
glet (right-handed) neutrinos vr. However, the number 3 is not sacred. On the one 
hand, the 2 non-zero light-neutrino mass differences required by experiment could be 
explained with just 2 heavy right-handed neutrinos. On the other hand, 3 non-zero 
masses could be explained with the participation of any number TV > 3 right-handed 
neutrinos. From the bottom-up perspective, there are no experimental constraints 
on the number of right-handed neutrinos. For example, since they are singlets of the 
Standard Model (SM) gauge group, their presence below the GUT scale would not 
perturb the unification of the SM gauge couplings. It is certainly possible, and even 
appears quite plausible, that such SM singlets do not follow the family pattern of the 
SM fermions. 

Specific examples of scenarios with many right-handed neutrinos are provided by 
string models, and some recent string constructions motivate seesaw models with up 
to 0(100) right-handed neutrinos [5], [6]. The reason is that string models contain 
abundant SM singlets that often have (non-renormalizable) couplings to the SM lep- 
ton doublets and large Majorana masses, which are the ingredients necessary for the 
seesaw mechanism. This is just one of many scenarios that motivates a generalization 
of the conventional seesaw to models with many right-handed neutrinos. 

In this Letter, we explore some of the phenomenological features of such scenarios. 
Specifically, we comment that the seesaw mass scale may rise to larger values than 
in minimal schemes, that the magnitudes of lepton-flavour- and CP-violating effects 
may increase, and that leptogenesis bounds on the reheating temperature and the 
lightest neutrino mass may be relaxed. 

2 Formulating the Seesaw with many Right- 
Handed Neutrinos 

The supersymmetric seesaw mechanism is described by the superpotential (see e.g. [7]) 

W = Yl^ eilj + Yi^ 0" Nitj + ^M.-feAT,- , (2.1) 

where (f)^''^ and £i {i = 1,2,3) are the Higgs and lepton doublets, Cj are the charged 
SU(2)-singlet leptons and Nj {1 < j < n) are some heavy Standard Model (SM) 
singlets with Majorana mass terms Mj^. The Yukawa couplings form an n x 3 matrix, 
while the Majorana mass terms form an n x n matrix. The resulting effective mass 
matrix for the left-handed neutrinos is given by 

Meff = -{vsinpfYjM-'Y^ ^ - ^11^ S, (2.2) 
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where v = 174 GeV, tan/3 is the ratio of the Higgs VEVs, is the "effective seesaw 
scale" and 5 is a 3 x 3 texture whose largest entry is of order one. 

As we now discuss, low-energy physics is quite sensitive to the number of right- 
handed neutrinos n participating in the seesaw mechanism. In particular, the effective 
seesaw mass scale may depend on n: 

oc n-"" , (2.3) 

where < x < 2 is a function of the texture. We consider the following limiting cases 
(in a particular basis defined, for example, by the Proggatt-Nielsen charges): 

(1) all entries of M and contribute coherently to Mes with similar magnitudes, 
so that oc 

(2) the magnitudes of the different entries of M, Y^, are similar, but the complex 
phases are randorD0, so that M* oc n~^; 

(3) only a limited number of right-handed neutrinos provide significant contribu- 
tions, so that oc nP. 

Typical realistic scenarios correspond to intermediate situations. For instance, 
if M is approximately diagonal with similar eigenvalues and all the Y^'' are similar 
in magnitude but have arbitrary phases, M^, oc n~^/^. As a specific example, we 
note that the stringy model of Ref. [5] corresponds to x between and 1. 

In most cases, the masses of the left-handed neutrinos grow with n. Therefore, 
in order to keep them at the same values as in the conventional seesaw with 3 right- 
handed neutrinos, one has either to increase the Majorana masses or to decrease the 
Yukawa couplings. These possibilities differ in their phenomenological implications. 

Consider, for simplicity, an MSSM scenario with universal soft supersymmetry- 
breaking scalar masses at the GUT scale, as motivated by minimal supergravity 
(mSUGRA), complemented by such a seesaw with n right-handed neutrinos. The 
renormalization-group (RG) running between the GUT scale and the scale at which 
the heavy neutrinos decouple that is due to the neutrino Yukawa couplings induces 
additional flavour-violating soft terms [8], 

where m? is the slepton mass-squared matrix. These corrections are given by [9]- [14] 
Sml^-^{3ml + Al)Y^LY,, 



^In this case, the dependence of on n is obtained using the "random walk" result that 

N 

where Ui is a random integer. 



Here = ln(MGUT/Afi) ^ij with Mj being the heavy neutrino mass eigenvalues. 

We first note that the dependence on the Majorana masses is only logarithmic. 
Thus, increasing the scale of the Majorana masses does not affect these corrections 
significantly. On the other hand, Y^Yy in general grows with the multiplicity of the 
heavy neutrinos, plausibly as a power law. Therefore, ignoring the logarithmic piece, 
we have 

5m\ , 5m\^ oc nP , (2.6) 



where the power < u < 1 depends on the texture. In the extreme case of n "coherent" 
contributions, y = 1 □; for a case with random phases, y = 1/2; finally, y = when 
only a few heavy neutrinos contribute significantly. Clearly, fiavour-violating effects 
are expected to be enhanced, in general. In particular, the branching ratio for radiative 
h — *■ Ijl decays [9], [11], 

r(/, ^ / oc Q^mf, tan2 (3 , (2.7) 

where rh characterizes the typical sparticle masses in the loop and a is the fine 
structure constant, is enhanced by n^^ which can be a large factor (up to 4 orders of 
magnitude) H. 

Alternatively, if instead of increasing the scale of the Majorana masses, one de- 
creases the Yukawa couplings, no enhancement of lepton-flavour-violating (LFV) is 
expected, in general. This is because M^q and 6m\ j^j^ both depend on Yi, quadrati- 
cally. 



3 Phenomenological Constraints 

The requirement of perturbativity restricts the magnitude of the loop corrections, 
which translates into 

This is a rather weak constraint and leaves open the possibility that some entries of 
YuLYy could be as large as 10. We note that there are no (very) large logarithms in 
the loop corrections since the seesaw scale is quite high. 



^ Unlike in Eq. l|2.3p . here we work in the basis where the Majorana mass matrix is diagonal. Thus, 
there are n contributions in the sum which we assume to be similar in magnitude in this basis. Note that 
the basis change may, in the case of many vr^s, change the orders of magnitude of the couplings. 

^We also note that, since Y^Y^ is larger than in the usual case, the presence of many right-handed 
neutrinos improves the MSSM gauge coupling unification at two loops [15]. 
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The main constraints are imposed by the lepton- flavour-violating branching ratios 
[16]- [18] 

BR(^ ^ 67) < 1.2 X 10"" , 
BR(r ^ 67) < 1.2 X lO"'^ , 

BR(r ^ ^7) < 6.8 (4.5) x 10"*^ . (3.2) 

Supersymmetric flavour-violating effects can be expressed in terms of the "mass in- 
sertions" 5lr = 5m\^/m? and Sll = 6m\/m?' [19], where m is the average slepton 
mass, such that [20] 

BR{k - 1,1) = |6 + 6 ^Ir^ + {L^R). (3.3) 

Here ^1^2 are functions of SUSY masses whose explicit form is given in Ref. [20]. For 
electroweak-scale sparticle masses, the resulting constraints on the mass insertions 
are [20]: 

l^iil < few X 10-3 , iJf^l < 10-6 , 

\5f^\<l^-^ , |5i3j|<10-2. (3.4) 

By means of Eq. (l2.5l) . these bounds are translated at low tan/3 into 

{Y}LY,)i2 < 0.1 , (3.5) 

with the other entries being essentially unconstrained. 

The interpretation of this constraint depends strongly on the magnitudes of the 
Yukawa couplings and the details of the texture. For smaller Yukawa couplings, the 
number of RH neutrinos can be very large. However, for couplings of order unity, 
only a few RH neutrinos are allowed, unless there are cancellations. In general, even 
though (yJlY^) 12 is rather small, other matrix elements can be significant and lead 
to observable effects. 



4 CP- Violating Phases and Electric Dipole Mo- 
ments 

The Yukawa couplings and the Majorana mass terms are in general complex, leading 
in general to CP violation. The number of physical CP-violating phases in the high- 
energy theory can be determined by parameter counting. Initially, Ye has 9 complex 
phases, Yi, has 3n phases, and M has n(n -|- l)/2 phases. On the other hand, the 
"flavour" rotation symmetry 

UlXUeX Un, (4.1) 
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acting as 

Ye ^ Ul Ye Ui , 
Y, ^ U^^ Y, Ul , 

M ^ Uj^MU*^, (4.2) 

allows one to eliminate some of these phases by field redefinitions. We note that Ui 
and Ue contain 6 phases each, while Un has n{n + l)/2 phases 0. Thus, this flavour 
rotation leaves 3(n — 1) physical phases, which we can identify explicitly in a specific 
basis. Consider the basis where 

= real diagonal , 
M = real diagonal , 

= arbitrary . (4.3) 

This basis is defined only up to a diagonal phase transformation 

Ui = Ue = dmg{e'^\e''^\e''^^) , (4.4) 

which acts on the neutrino Yukawa couplings as Y'J-' Yy^ e^'^'^ . The physical phases 
are invariant under this residual symmetry and are given by 

Arg(y;^ yj=^ *) , (4.5) 

(where no summation over j should be understood). Clearly, there are exactly 3(n — 1) 
such phases, and each one may play the role of a "Jarlskog invariant" (see [21], [22] 
for a discussion). 

Complex phases in the neutrino Yukawa couplings induce CP-violating phases in 
the soft terms due to the RG running, which in turn contribute to lepton electric 
dipole moments (EDMs) at low energies [23]- [27]. The relevant flavour objects are 
the 3x3 matrices yJy^ and yJlY^, such that [23], [27] 



dj oc 



Y^Yy , Yy LYy 



(4.6) 



where di is the EDM of the i-th charged lepton. This expression is proportional to 

(no summation over i) which makes it clear that only 
the reparametrization-invariant phases (|4.5I) are involved. Since the summation over 
n RH neutrinos appears twice, the EDMs grow as 

di cx . (4.7) 

Thus, they may be enhanced by several orders of magnitude in the multi-neutrino 
case. 



"'In the pure Dirac case, one of these phases is irrelevant since an overall phase redefinition leaves all 
flavour objects invariant. 
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At large tan/3, a different flavour structure appears [26], [27], namely: 

di oc im (yj/y, m2 yj^y,) , (4.8) 

where f,g are diagonal matrices depending on the Majorana masses and nii are the 
charged lepton masses. The conclusion, however, remains the same, and the EDMs 
grow with the multipUcity of the states as di oc n^^. 

The order of magnitude of the induced EDMs can be estimated using [27]: 




where we have taken ln(Mj/Mj) = 0{1) and CP-violating phases that are of order 
unity. The current experimental limit on the electron EDM of 10~^^ e cm is saturated 
for (yjy^) 11 ~ 10 if the superpartners of SM particles have masses at the electroweak 
scale. For large tan/?, the above expression acquires an additional factor (tan/3/10)^ 
which leads to a further enhancement of the EDMs. 

It is important to remember that these estimates are very sensitive to other CP- 
violating phases. For example, if the phase of the //-term is as small as IQ-^-lQ-^, it 
will dominate the SUSY contributions to the EDMs (for a recent discussion, see [28]). 



5 Numerical Example 



In this Section, we illustrate the multi-neutrino scenario with a numerical example. 
We suppose that, in the basis where the charged lepton mass matrix is diagonal and 
positive, the mass matrix for the light neutrinos is given by 

/ 0.003 0.003 0.003 \ 
Meff = 0.003 0.028 -0.022 

V 0.003 -0.022 0.028 / 

in eV units. This mass matrix has eigenvalues (0.05,0.01,0) eV and is diagonalized 
by a tri-bimaximal [29] PMNS transformation. In our convention, the (1,1) entry of 
Meff corresponds to the r-neutrino, (2,2) - to the muon neutrino and so on. 

We consider two seesaw realizations of this light-neutrino mass matrix: one with 
10 right-handed neutrinos and another with 2 right-handed neutrinos. The corre- 
sponding heavy-neutrino mass matrices are generated randomly under the conditions 
that the Yukawa couplings be of order one and that most of the entries of M^^ be 
similar in magnitude, while reproducing the correct M^s- For simplicity, we choose 
real matrices. 

In the case of 10 right-handed neutrinos, the inverse Majorana mass matrix M^^ 
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is given by 
/2.I 



-1.58 
1.58 0.13 



0.02 
0. 

0.02 
0.92 
0.15 
3.17 
0.22 
1.65 



5.09 

0.14 

-0.33 

3.51 

-0.16 

-0.35 

-3.41 

1.61 



0.02 

5.09 

0.83 

-0.09 

2.85 

-0.97 

-0.91 

2.4 

-0.25 
4.45 



0. 

0.14 

-0.09 

-4.19 

-3.79 

1.11 

-1.03 

0.45 

-1.5 

2.79 



0.02 0.92 
-0.33 3.51 



2.85 

-3.79 

-4.6 



-0.97 

1.11 

-0.12 



0.15 

-0.16 

-0.91 



3.17 

-0.35 

2.4 



0.22 
-3.41 



1.65 
1.61 



-0.12 0.25 

-4.42 0.35 

1.91 0.41 

-0.45 0.74 

-0.31 2.58 



-1.03 0.45 

-4.42 1.91 

0.35 0.41 

-3.07 2.61 



-0.25 4.45 
-1.5 2.79 



2.61 
2.24 
-2.64 



1.9 

4.51 

-0.15 



-0.45 

0.74 

2.24 

4.51 

4.91 

0.33 



-0.31 

2.58 

-2.64 

-0.15 

0.33 

2.15 



in units of —3.3 x 10 GeV ^, where we have assumed sin/3 1. The corresponding 
Yukawa couphngs YJ" are 

0.51 0.06 -0.32 0.32 -0.83 -0.36 1.23 0.41 -0.35 0.48 
0.37 0.1 1.36 0.62 -0.82 -0.66 0.06 -0.02 -0.78 0.11 
-0.09 0.89 0.13 -0.55 0.08 -0.25 -0.04 0.68 0.11 -1.1 
The Majorana mass eigenvalues are Mi =(2.6, 2.7, 3.1, 4.0, 5.0, 8.1, 9.3, 10.5, 18, 
58)xl0^^ GeV with the geometric average scale 7.2x10^^ GeV. In the basis where the 
Majorana mass matrix is diagonal, the Yukawa couplings remain of order one. 

In the case of 2 right-handed neutrinos, we nevertheless use a 3x3 matrix notation, 
keeping in mind that has rank two, so that only two right-handed neutrinos 

contribute. We take 



M 



-1 



/ 11.28 -8.82 21.86 

-8.82 6.91 -17.21 

\ 21.86 -17.21 44.39 

in units of -3.3 x lO^^^ GeV-\ and 



1/ 




(0.5, 56)xlO^'' GeV with the geometric 



/ -0.02 

0.83 
V 0.54 

The Majorana mass eigenvalues are Mi 
average scale 5.2x10^^ GeV. 

As expected, the Majorana mass scale is somewhat higher in the first case, 
while the Yukawa couplings are similar in magnitude. For the comparison of their 
lepton-flavour- and CP-violating effects, the relevant quantity is YuLY^, where 
L = diag[ln(MGUT/-^j)]- Ignoring the RG running of Y^ [30], we have 
/ 7.47 -0.39 -1.65 \ 



V 



-0.39 9.77 -0.58 
-1.65 -0.58 8.4 



^More precise numbers are available from the authors. 
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for the 10 right-handed neutrinos case and 

/ 0.67 1.19 0.15 \ 
Y^LY^= 1.19 2.37 0.01 

V 0.15 0.01 0.29 / 

if there are 2 right-handed neutrinos. Here, as before, the (1,1) entry corresponds to 
the T-T matrix element and so on. 

We find that not all of the entries are enhanced by the same factor in the 10 
right-handed neutrino case, but the phenomenological consequences are dramatic! 
For example, the r — > 67 branching ratio increases by a factor of 100 compared to 
the 2 right-handed neutrino case, while that for /x — > 67 is enhanced by a factor of 10'^. 
This result is of course specific to the above texture, however it is generally true that 
the lepton-flavour-violating transitions are more significant in the multi-neutrino 
case. Although here we have considered a CP-conserving Ansatz, similar statements 
apply in general to the EDMs as well. Since the magnitudes of the diagonal entries 
increase by factors of 10 or so, the insertion of CP phases of order unity would result 
in increases of EDMs by two orders of magnitude. 

The above example illustrates the main features of the multi-neutrino scenario. 
This particular, the 10 right-handed v texture would be consistent with the BR(/i 
67) constraint only for slepton masses of ~ 500 GeV or more. The corresponding 
predictions for BR(r — > 67) and BR(r — > fi'^) are at the 10"^'^ — 10^^^ level. For 
low tan/3 and CP-violating phases of order unity, the electron EDM is expected to 
be of order 10~^^ e-cm and the muon EDM a factor of 10^ larger. Thus, the most 
promising observables for this texture would be BR(/i — > 67) and de (see [31] for a 
discussion of the experimental prospects). These flavour-dependent effects would be 
accentuated further for more right-handed neutrinos. 



6 Comments on Leptogenesis 

A number of features of leptogenesis [32] with many right-handed neutrinos have 
been studied by Eisele in Ref. [33]. Here we discuss only a few of the most important 
differences of this scenario compared with the standard scheme [§. 

The out-of-equilibrium decay of the lightest right-handed neutrino creates a lep- 
ton asymmetry 

VL = e 6 , (6.1) 

where n and s are the number and entropy densities, respectively; 6 is the washout 
parameter characterizing the fraction of the lepton asymmetry surviving the washout 



'Recent reviews of leptogenesis can be found in Ref. [34]. 
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effects, and e is the CP-violating asymmetry: 

r -r 



r + r ' 



(6.2) 



with r being the decay rate of the singlet right-handed neutrino into light leptons 
plus the Higgs boson, and F being the decay rate for the CP-conjugate process. In 
the case of hierarchical Majorana masses, the CP asymmetry is dominated by the 
decays of the lightest singlet, and 

3 Ml im [y.Mj^yJ]^^ 

where Mi is the lightest Majorana mass. Davidson and Ibarra have derived an upper 
bound on ei in the 3x3 case [35]. Here we consider the corresponding bound in the 
multi-neutrino scenario. Working in the basis where Meff is diagonal and real, we 
first note that the asymmetry is maximized for Y^^ = \\Yii'\^. Then, the relevant ratio 
satisfies 

T.^m,\Yu\'' 



|2 



< "i3 , (6.4) 



where is the mass of the heaviest left-handed neutrino. The maximum is achieved 
for yii,yi2 = 0. Thus, we obtain a generalized bound for an arbitrary number of 
right-handed neutrinos: 

3 Ml 
8^ 



l^il ^ ^ 7W ■ (6.5) 



This is identical to the "weaker" version of the Davidson-Ibarra bound in [35], which 
was also found in Ref. [36]. For the 3x3 case, a stronger bound is valid, which is 
obtained from the above expression by replacing with m^ — mi. This stronger form 
does not apply to the n x n case, and one can construct examples with yii,i2 = 0, 
Arg(yi3) = 7r/4 which saturate (|6.5p . Technically, this happens because the Casas- 
Ibarra parametrization [11] involving orthogonal matrices does not apply to n x 3 
Yukawa couplingil|. 

Thus, only the weaker version of the Davidson-Ibarra bound holds in the general 
case. Needless to say, this distinction is only relevant for a degenerate light-neutrino 
spectrum. In the 3x3 case, the bound on e tightens as the overall scale of light neutrino 
masses increases because m^ — mi ~ Am^tm/2"i3- For many right-handed neutrinos, 
the only relevant bound is (16. 5p which relaxes as ms increases. This inequality can 
be interpreted as a lower bound on Mi, which in turn implies a lower bound on the 
reheating temperature Tj-eh [35]. If we take ms ~ 0.5 eV, the Davidson-Ibarra bound 



on Tj-eh relaxes by an order of magnitude 

Treh > 10^ - 10^ GeV . (6.6) 



^The n x 3 matrix R entering the Casas-Ibarra parametrization satisfies if^ R — 1 but not RR^ — 1. 
^For systematic studies of bounds on T^eh, see [37], [38]. 
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This ameliorates somewhat the problems of gravitino over-production [39] and moduli 
destabilization at high temperature [40]. 

Another interesting consequence of the multi-i^/j scenario is the relaxation of the 
cosmological bound on the lightest neutrino mass mi. The requirement of out-of- 
equilibrium decay of the lightest right-handed neutrino [42], Fi < H\t~m-^, implies 

mi = (y.yj)ii-^ < 5 X 10-=^ eV . (6.7) 

In the 3x3 case, mi < mi (see, for example, [43], [35]). In the multi-i^ij case, this is 
no longer true, since mi receives n contributions, mi oc {YjY,y/Mi)ii, and one can 
easily have mi > m-i. Thus, the constraint on the lightest neutrino mass is relaxed. 

We see that the presence of more than 3 right-handed neutrinos disrupts the usual 
relations between the light neutrino masses and leptogenesis. In the above arguments, 
we have made the usual assumption that the lightest i'r dominates in leptogenesis. 
We have assumed that the extra vrs play a passive role, in the sense that they do 
not contribute to leptogenesis, yet affect the light neutrino masses. 

Further possible effects of many right-handed neutrinos have been discussed in 
Ref. [33]. These include, for example, the possibility that more urs are responsible 
for leptogenesis. Such effects can be studied on a model-by-model basis. 

7 Conclusions 

We have explored some of the phenomenological features of the seesaw mechanism 
with many right-handed neutrinos. Such a generalization is allowed by the absence 
of experimental constraints on the number of heavy right-handed neutrinos, and is 
motivated by string constructions. We find that, in such models, the Majorana mass 
scale can be significantly higher than that in the traditional scheme with 3 right- 
handed neutrinos, and have demonstrated with a specific example that the magni- 
tudes of lepton-flavour- and CP-violating transitions can be enhanced by orders of 
magnitude. We also find that certain constraints on leptogenesis are sensitive to the 
number of and can be relaxed. This applies, in particular, to the bounds on the 
reheating temperature and the lightest neutrino mass. 
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